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Abstract 

We discuss several aspects of the cosmological evolution of moduli fields in het- 
erotic string/M-theory scenarios. In particular we study the equations of motion 
of both the dilaton and overall modulus of these theories in the presence of an ex- 
panding Universe and under different assumptions. First we analyse the impact 
of their couplings to matter fields, which turns out to be negligible in the string 
and M-theory scenarios. Then we examine in detail the possibility of scaling in 
M-theory, i.e. how the moduli would evolve naturally to their minima instead 
of rolling past them in the presence of a dominating background. In this case 
we find interesting and positive results, and we compare them to the analogous 
situation in the heterotic string. 
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1 Introduction 



For a number of years now, our knowledge of the string world has greatly improved with 
the discovery of M-theory as the origin of all the known perturbative string theories 0. 
In particular, the D — 11 supergravity (SUGRA) limit of this theory (also known 
as M-theory itself) has been extensively studied and, by now, we have quite a deep 
understanding of many aspects of it. 

However the issue of how to stabilize the moduli in these theories, and the role 
they play in cosmology, is still an open one. This is because, in order to generate a 
potential for the moduli, we must rely on non-perturbative physics which is not under 
our control. Still, a few general results can be obtained, which encourage further inves- 
tigation. In particular, supersymmetry (SUSY) breaking by gaugino condensation 
in the D = 4 effective theory works in a similar way as in the old heterotic string ||. 
In fact it is possible to build many scenarios where multiple gaugino condensates con- 
tribute to create a scalar potential for the moduli S and T that has minima at the 
phenomenologically desired values (Re S, Re T ~ O(20) in Mp units) [|J. Also the 
fact that the gaugino condensates depend exponentially on the gauge kinetic functions 
(~ e~ a °-f a , where a labels each condensate, a a are related to the 1-loop beta function 
coefficients and f a are the gauge kinetic functions), with 47r/ a = S — (3 a T 1 instead of 
fa = k a S as it was the case in the heterotic string, gives rise to the appearance of flat 
directions in the scalar potential, opening up the possibility of inflation along those. 
This was thoroughly studied in Ref. ||. 

Despite of these promising results, there is a second issue we want to address in this 
letter, which is the dynamical evolution of the other moduli, i.e. those which do not 
have a flat direction and therefore would suffer from a 'runaway problem', given the 
exponential nature of their potential. This problem, first pointed out by Brustein and 
Steinhardt || in the context of the heterotic string, has been studied along the years 
and several mechanisms have been proposed in order to alleviate it (see Refs 0^9|). 
Our goal will be, exploiting the similarity between the heterotic string and M-theory 
scenarios pointed out above, to examine this mechanism for dynamically stabilizing 
moduli in the context of heterotic M-theory. 

In order to do so, in section 2 we revise the status of the runaway problem in the 
heterotic string, with particular emphasis on the recently proposed mechanism of Huey 
et al. ||, and its applicability to string models. In section 3 we discuss the M-theory 
setup: first we study the existence of scaling solutions in the presence of a dominating 
background and then we consider the above-mentioned couplings of the moduli to 
matter. Finally, we conclude in section 4. 

2 Scaling solutions revisited 

Let us briefly review the issue of scaling solutions in gaugino condensation models. It 
is by now a standard result in cosmology that scalar fields with an exponential scalar 
potential have scaling solutions that are global attractors [|I~0| -|T2|| . In these scaling 
solutions the field will evolve with the same equation of state as the background energy 
density. As was explained in Ref. M , the dilaton, S, in the heterotic string theory is a 
perfect example of such scaling behaviour. This scalar particle has an exponential type 
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potential which is steep enough to make the field roll past its minimum. However, if the 
energy density of the Universe includes some barotropic fluid coupled to the dilaton 
only gravitationally, this field will reach a constant velocity at late times, and will 
evolve slowly down the exponential slope. This behaviour, however, can only stabilize 
the dilaton if it starts its evolution to the left of its minimum. 

More recently a new mechanism has been proposed || which could help stabilizing 
moduli with exponential potentials both in their strong and weak coupling regimes. It 
relies on the existence of couplings of these moduli (generically denoted by 0, with an 
exponential type potential Vo (</»)) to matter fields (which we will call C). These can 
occur either through the kinetic terms (which, following Ref. |§, we parametrise as 
f(4>)d^Cd^C) or through the scalar potential (Vi ~ g(<f))C n ). Assuming the C field to 
be homogeneous, we define its energy density and pressure to be respectively, 

Pc = fC 2 + V, , (1) 
Pc = fC'-V, . (2) 

The equations of motion for cf) and p are then given by 

+ 3iJ0-^ + — = 0, (3) 
Pc = -(^H+ (p c + Pc ) + 7^(Pc- Pc) , (4) 

with H 2 = |((/> 2 /2 + Vq + pc)- To be precise, we will consider the case where the C 
field has coherent oscillations about its minimum. In this case, the energy density and 
pressure, averaged over one oscillation, will obey the equation of state pc = wcpc, 
with wc = (to — 2)/(n + 2) from the virial theorem. We do this in order to have a 
simple barotropic fluid as a background. A similar situation could be obtained if we 
considered instead the field C to be in thermal equilibrium 0. 

With these assumptions, the pressure derivative in Eq. (|3|) becomes 

dp c (<P,C,C) ( ldf ,ldg\p c 

and Eq. fl4|) can be exactly solved to give 

where a is the scale factor of the Universe and the subscript o denotes initial values of 
the corresponding functions. Altogether we see that the new term —dpc/dip can be 
interpreted as a contribution to the total scalar potential of the field. In other words, 
the coupled system will be equivalent to a single field system with a total potential 
given by 

Kot = W) + K ff (0) , (7) 

with 

Vb(0) ~ e"^ (8) 
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Figure 1: 

(a) Plot of the total scalar potential, given by Eq. (j^), as a function of the dilaton cj) (in Mp units) in 
a two condensate model with groups SU(6)2 x SU(7)g and Tr = 1.23 (see footnote 2). The different 
curves correspond to different times -or number of e-folds N- in the evolution of the system, (b) 
Evolution of the dilaton (<ft), as in Eq. (|ll|), as a function of N = ln(a), for the same example as in 
Fig. la and initial conditions given by 0o = 13 (dashed), 4>q = 4 (solid). 



and 

VM) = p c (0) • (9) 

Notice that pc is now given by Eq. (|6]), an explicit function of time (or the scale factor, 
a), as opposed to the pc we had before integrating out the C field. They have both the 
same value of course, but their partial derivatives with respect to (p will be different. 
In particular, it is easy to show that 

dpc{4>,a) dp c ((f>, C, C) 

the result we have used to obtain Eq. ([|). This effective potential is represented in 
Fig. la, where we show the shape of Vtot as a function of <ft, for several values of 
N = ln(a) with wq = 1/3, for the particular case of /(</>) = g(4>) = 1/0 considered 
in Ref. 0. Here V is given by a two condensate model analogous to Eq. (4) of Ref. 0. 
We can see that, at earlier times (i.e. for N small), the potential has a growing slope 
to the right of the minimum, whose position changes with N. At late enough times 
(i.e. for N large) this effective contribution fades away and we recover the original 
exponential potential Vq(4>). As we see in Fig. lb, even for initial conditions in the 
very weak coupling regime (characterised by large values of 0), the field rolls towards 
its minimumF]. 

Despite of how promising this mechanism looks as a way of stabilizing moduli, it is 
very sensitive to the particular forms of the functions / and g (recall / = g = 1/0 in 
the above example). Therefore it is a natural step to proceed with the study of string 

1 It can be argued that, in the weak coupling regime, it would be inconsistent to consider the C field 
in thermal equilibrium with the background. However this is beyond the scope of this paper which is 
to show that in a string motivated scenario the effect of matter fields is negligible in any case. 
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inspired models and, in particular, we shall devote some time to consider the dilaton 
of the heterotic string and its couplings. 

The first thing one must take into account is that, in the context of the heterotic 
string, the dilaton S is not a canonically normalised field, i.e. its kinetic term is given 
by Kggd^S d^S, with K s § = d 2 K/dS3S and the tree-level Kahler potential given by 
K = — ha.(S + S). Therefore its equation of motion reads 

where K s8s = dK s§ /dS and H 2 = l(K s§ S 2 + K c0 C 2 + V(S,C)). It must be said, 
however, that the term proportional to S 2 is negligible in the region of values of S 
studied here, and therefore this equation is very similar to the simpler Eq. (^|) described 
previously. The couplings to matter, contained in dpc/dS, are determined by the 
structure of the N=l SUGRA Lagrangian that describes this string-inspired modelQ. 



The kinetic term for C reads QT3 1 



K^CPC = ( 1 + -|r ) d.O < < . ( 12) 



2T R V 6S R 

where e = <5gs/47t 2 is the so-called Green-Schwarz coefficient, which is generally small, 
and S R = ReS, T R = ReT. Therefore in the parametrisation introduced before Eq. ([!]), 
f(S R ) is given by 
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e 



Finally, the form of the scalar potential is given by 

V = e K i^(Wi + KiW)(Ki)~\w j + K j W) -3\W\ 2 ^ , (14) 

where sub(super)scripts denote derivatives of W and K with respect to S (S) and 
C (C), and W ~ J2a e ~ aaS + C 3 is the superpotential coming from multiple gaugino 
condensation. The couplings between the two fields S and C arise from the term 
e K \KcW + Wc\ 2 / K C (j. For C small, Wc ~ C 2 will dominate and we can identify 

g(S R ) = e K /K cC = — 1— . (15) 

Therefore this coupling is of the form 1/S R . 

We can now solve Eq. (JTT]) , and the result is that S R has the same runaway behaviour 
that it would have had in the absence of the couplings to matter. Let us briefly sketch 
why the mechanism of Ref. does not work here: consider the general form for the 
couplings / = a + b/S R , g = c + d/S R , with a, b, c, d real and positive. A slope to the 
right of the minimum, such as those in Fig. la, will appear if, in Eq. (O), we require 



2 Note that, throughout this section in which we focus on the evolution of the dilaton, we are 
assuming that the T modulus, whose vacuum expectation value (VEV) represents the overall size of 
the compactihed space, is fixed to its value at the minimum. We will include the T modulus explicitly 
when considering heterotic M-theory scenarios in the next section. 
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—dpc/dS > to balance dV/dS which is invariably negative (recall that V ~ e~ aS ). 
From Eq. (^) is it easy to translate this condition, for positive Wc, into either of these 
two 

i) ad(l — wc) — bc(l + wc) < 

(16) 

2bdwc 



ii) ad(l — wc) — bc(l + wc) > and Sr. < 



ad(l — wc) — bc(l + wc) 



In Ref. 0, indeed wc = 1/3 is positive and / and g are such that a = c = 0, that 
is this choice fulfils condition i) and their field is always trapped at the minimum. 
In the heterotic string, however, c = and a ^ 0, so we are in case ii), which implies 
that the dilaton will be trapped only for initial conditions 5r < i/a = e/6 -C 1, for any 
values of wc ^ 1. This corresponds to the strong coupling regime and covers a very 
small fraction of the parameter space. For completeness we can also consider the case 
in which — 1 < wc < 0, and now we will have only one condition for —dpc/dS > 0, 
which is 

ad(l - w c ) - 6c(l + w c ) < and S R > — r . (17) 

ad{± — wc) — bc[l + wc) 

It is easy to see that, again, for the heterotic string where c = 0, this condition can 
never be fulfilled. 

Therefore to end this section we should conclude that, concerning the probability of 
stabilizing the moduli in the Early Universe, the estimates of Home and Moore [|7| still 
represent the most optimistic situation. They argued that the motion of the moduli 
follows a chaotic trajectory and, out of the finite volume of possible initial conditions, 
12-14% of it corresponds to cases where the dilaton will end up at its minimum. How- 
ever, in Ref. |TJj] it was pointed out how, in their scenario, the inhomogeneous modes 



would soon come to dominate the energy density of moduli fields. Finally, if we con- 
sider the presence of a dominating background, then the estimates of Ref. M remain 
unaltered in the presence of above-mentioned couplings to matter fields. 



3 Dynamics of the moduli in M-theory 

In this section we will try to apply the same ideas, namely the evolution of moduli fields 
in the presence of some background matter, to study the evolution of the real parts of 
the S and T moduli in heterotic M-theory. From a phenomenological point of view, one 
of the main differences between the two models is that the gauge kinetic functions in the 
hidden wall, f a (whose VEV determine the gauge coupling constants) are now a linear 
combination of the string dilaton and modulus fields, i.e f a = (S — /3 a T)/(47r) where j3 a 
are constants that depend on the details of the model and are usually of order one |Tj| 
(in particular, in all the examples shown in this paper we have fixed (3 a = (3 = 1/2). 
Assuming once again that gaugino condensation is the source of SUSY breaking, the 
scalar potential will have an exponential like profile, as in the string case, along the 
direction defined by $_ = S — /3T. Along its orthogonal direction, $ + = (3S + T, 
the potential will be almost flat, the only dependence upon this variable coming from 
the Kahler potential K. We will consider potentials where a minimum is generated 
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along the $_ direction with two gaugino condensates. This defines a 'valley' in the <3> + 
direction with an inverse power law profile, and the global minimum is obtained with 
a third condensate, or a non-perturbative correction to the Kahler potential @,[5|]. 

Let us start with the simplest case, when the S and T fields have no direct coupling 
to the background fluid. As we have just said we have essentially two major directions 
in our scalar potential, an exponential slope for the $_ field and an inverse power-law 
for the $ + field. In the exponential part of the potential we expect the field $_ to 
dominate the evolution, behaving in a similar way to the dilaton in the heterotic string 
case ||. However, in the valley defined by $_ = $ mm , we expect the evolution to be 
driven by the $ + field along the valley. Therefore the stabilization of the fields will 
rely on a two-step process, first the evolution into the valley, and then the stabilization 
along the valley to the minimum (see Fig. |^). 

We will analyse the behaviour of the fields in this model using a numerical simulation 
with a specific example. The equations of motion of the two fields, with no direct 
couplings to the background, is given by 

S + 3H S + ^ils> + J-™ = (18) 

f + 3Ht + ^ + J_ev = , (19) 

with the Hubble constant now being given by 

H 2 = ^(K s§ SS + K Tf ff + V + p B ) , (20) 

where pe is the background energy density. We evolve both fields simultaneously, since 
the flat directions in the scalar potential do not coincide with the ones in the Kahler 
potential. As detailed in Ref. ||, the scalar potential is explicitly given by 

V = e K \-^\Ws + KsW\ 2 + -^\W T + K T W\ 2 -3\W\ 2 ) , (21) 
[K ss K TT J 

where W, the superpotential for multiple gaugino condensation, reads as 

W = J2C a e' a ^-, (22) 

a 

with C a and a a being constants related to the one-loop beta-function coefficients of 
the corresponding condensing groups. We will consider for the numerical examples the 
case of two condensates, SU(3) x SU(4), with 8 pairs of matter fields transforming as 
(4,4). The Kahler potential is given by 

K(S R , T R ) = - logOS + S) - 3 log(T + f ) + K np (S R ) , (23) 

where the first two terms are the tree level result and K np accounts for non-perturbative 
corrections. We will use for the numerical examples the ansatz for K np introduced in 



Ref. 17 



K np = -£= log (l + e-^-VSo)) (24) 
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Figure 2: 



Contour plot of the scalar potential V in terms of Sr, Tr. The thick lines represent the trajectories 
of the fields in the presence of the corresponding backgrounds, in particular the upper right line 
corresponds to no background present. The parameters used for K np are B — 120 and D = 0.13. One 
can see that when a background fluid is present, the range of initial conditions allowing stabilization 
is highly improved. 

where we fix So = 19, and for each particular value of B we choose D so that the 
cosmological constant is equal to zero. 

Let us then start analysing the evolution of the fields outside the valley in the scalar 
potential. In Fig. we show the numerical evolution of the fields' equation of state, 
w<s>, as a function of before the fields reach the valley at $™ m = 11.07. A set of 
different initial conditions and different types of background fluid were considered. 

It is easy to see from this plot that the behaviour of the fields outside the valley is 
very similar to the case of the dilaton in the heterotic string case. That is, they will 
quickly reach a scaling solution in which their energy density scales with the background 
with roughly the same equation of state (i.e. w$ ~ wb = Pb/pb)- As in the string 
case, this will be enough to slow down the evolution of the fields and allow them to 
settle in the valley for a wide range of initial conditions. However, note in Fig. ^| that 
the evolution of the fields is not orthogonal to the direction of the valley. In other 
words, $_ and $ + scale simultaneously. The direction of the evolution is determined 
from the scalar potential gradient and the Kahler potentials of the two fields, and it is 
not possible to estimate it analytically. 

The main differences with the heterotic string case arise once the fields have settled 
in the valley. In order to be stabilized, they still have to evolve inside the valley 
towards the minimum. The valley corresponds approximately to a constant which 
we denote as $™ m , so that the evolution at this stage will be dominated by the $ + 
field. The effective scalar potential for the $+ field is power-law like. To be precise, 
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Figure 3: 

Equation of state of the system of fields $ + as a function of $_ before reaching the minimum 
at $™ m Rj 11.07. From top to bottom the different curves correspond to the following choice of 
background (%): 1/3, 0, -1/3, -2/3. 



the dominating term is the e K factor in Eq. QZTD, that is 
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+ $ min )($ + - /3$™ in ) 3 



(25) 



where Vq is a constant that can be estimated analytically (exactly for the D=0 case) and 
$ min is understood to be held constant. In this approximation, where $_ is constant, 
one can also write an effective kinetic term for the $ + field given in terms of the 
kinetic terms for the S and T fields. Namely, since S = + $_)/(l + (3 2 ) and 

T = ($+ - /?$-)/(! + (3 2 ), we have that 



/PKss + K, 



TT 



+ P 



2\2 



$ 2 



(26) 



The final result is a field which has a power-law type potential and a kinetic term that 
is inversely proportional to $+. If we transform it into a real field, a, with canonical 
kinetic terms, |d 2 , it is easy to check that a = V2 log$ + and therefore its effective 
scalar potential will be an exponential. Hence, we expect the behaviour at late times 
to be the usual scaling solution associated with exponential potentials. To be more 
precise, for large values of the field the effective potential, Eq. fl25D , goes as 1/$+, and 
therefore the effective scalar potential for o is V(a) = e" 2 ^. Without a minimum, 
that is for K np = 0, this is exactly the scaling solution one obtains numerically for late 
times. 

This is a nice result in itself, showing that the fields tend to join scaling solutions 
(albeit ones of different type) when in the presence of a background, both in the valley 
and outside it. However, the time the fields take to reach the scaling solution in the 
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valley is considerably longer. This means that they will usually reach the minimum 
of the potential before acquiring a scaling behaviour. This in itself does not stop the 
fields from being trapped at the minimum, since the background is already slowing them 
down considerably, but it does mean that an analytical approach is much more difficult, 
since the field reaches the minimum at an intermediate regime, between prescaling and 
scaling. 

We can have an idea of the range of initial conditions that stabilize the fields by 
looking at Table 1. First of all, if the value of p° B is too small, the fields will begin 
to evolve as if there was no background (that is, they will not stay in the minimum 
unless they start really close to it). Once the background becomes important, the 
actual value of p° B does not really affect the evolution, since the fields will just stay 
frozen until the energy densities of the fields and background become comparable. 
Therefore we chose to fix p° B = 10 20 in the examples given in Table 1. The value of wb, 
on the other hand, affects drastically the results. The lower it is, the easier it is for 
the fields to be stabilized. This is not too surprising, since a smaller wb corresponds 
to a universe where the potential terms dominate over the kinetic terms (and recall 
that the fields will mimic the background once they reach their scaling solutions). 
Moreover, different values of the (B, D) parameters in the non-perturbative Kahler 
will give slightly different answers. In short, the smaller the B parameter, the less fine 
tuned is the minimum in the scalar potential, and the easier it is to stabilize the fields. 
Finally, one can check that the further the fields are from the minimum when they hit 
the valley, the more difficult the stabilization becomes. The full evolution for specific 
examples where stabilization is achieved can be seen in Fig. |[ 
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Table 1: Examples of points in parameter space for which the fields are -or not- stabi- 
lized at the minimum. In all of them p° B = 10 20 and = 14. 



Let us finally turn to the case where the dilaton and modulus fields have direct 
couplings to the background fields. As for the string case considered in Section 2, we 
will again model the background as a scalar field C with a coupling /(Sr., Ir) in the 
kinetic terms, and g(S^, Tr) in the scalar potential. Its energy density and pressure is 
again given by Eqs (0,|j) and the dilaton and modulus equations of motion by Eq. ( |T8"D 
and Eq. (|T9|) respectively. The couplings are now functions of both fields |15| , namely 



/("Sr.j Ir) = ^ + ^- , (27) 
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The situation here is similar to the heterotic string case if slightly more involved. Again, 
dV/dS is always negative (except near the valley where the condensates cancel each 
other), and the S field tends to go to infinity. If the interaction term in Eq. (|TI|) is to 
cancel this natural tendency, then we require dpc/dS < 0. It is easy to check that, for 
positive S and T, this is equivalent to S < 2/3/3T. This only happens for S smaller 
than its minimum value, so if the initial conditions are larger than this, the field will 
always run away to infinity. For the T field, however, things do change. The derivative 
term dV/dT is always positive (again, except near the valley), and therefore the field 
will tend to go to zero (it will never become negative since its kinetic term diverges 
at zero). It is possible to check that dpc/dT has the same sign as pc and therefore, 
for positive background pressure, the coupling term balances the potential term in the 
equation of motion for T. Indeed, we have checked numerically that the T field will 
tend to go to a non-zero constant in the presence of these coupling terms. This can 
also be understood analytically if one looks for the asymptotic solutions in this regime. 

Finally let us add that this analysis was performed in the absence of non-perturbative 
corrections to the Kahler potential, K np (SR) (notice the form of Eq. (p8j) ). In the pres- 
ence of such corrections, that would only affect the behaviour of the S field, we would 
have to replace Eq. (p8|) by e Knp g(S^, Tr), and its derivative with respect to Sr would 
be now e Knp (dg /dS^ + gdK np /dS^). Given that dK np /dSR is also negative, this would 
not alter the conclusions we just reached about the evolution of S. 

In short, adding the interaction terms only improves the situation for the T field. 
In practice, this means that the fields still have to start to the left of the valley in the 
scalar potential if they are to become trapped at the minimum. We did not find any 
improvement in the allowed regions of parameter space presented in Table 1 when we 
included couplings with the matter fields. 

4 Conclusions 

In this letter we have addressed the dynamical evolution of moduli fields in several 
string/M-theory scenarios. In all of them the dynamics are provided by assuming 
gaugino condensation in the hidden sector of the theory as the source of SUSY breaking 
and, therefore, of a non flat potential for these fields. Following our previous work 
on the existence of scaling solutions for the dilaton evolution in the heterotic string, 
which provided a solution to the so-call 'runaway problem' pointed out by Brustein 
and Steinhardt, we investigated here the impact of a recent proposal to also stabilize 
moduli of considering the couplings of moduli to matter fields. After a thorough study 
of the issue in the context of the heterotic string we concluded that these couplings do 
not seem to affect stabilization for realistic string settings. 

Next we considered M-theory scenarios, where there is also a runaway problem 
associated, in this case, to the combination of dilaton and modulus, which has an 
exponential-type potential analogous to that of the dilaton in the heterotic string. Here 
we study the evolution equations for both $ + and $_ in the presence of a dominating 
background, in order to see whether we can also reach a scaling regime that will make 
the fields settle at their minima instead of rolling past them. It turns out that such 
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a regime, i.e. scaling, is achieved but at late stages of the moduli evolution, with the 
result that these fields are not stabilized while scaling, but in an intermediate regime. 
This makes it extremely difficult to perform an analytic study of the problem but, 
nevertheless, we have been able to determine a wide region of the parameter space 
(essentially defined by the initial positions of the fields, the type of background and 
the characteristics of the non-perturbative Kahler potential) for which the stabilization 
is successful. Finally, we considered the couplings of moduli to matter fields, and we 
concluded that the situation is the same as in the heterotic string case, i.e. these 
couplings do not contribute at all to improve the stabilization problem of the moduli. 
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